//o 7c 

//!/ - j?7 

i 


NASA AVSCOM 

Technical Memorandum 100180 Technical Memorandum 87-C-2 


On Dynamic Loads in Parallel 
Shaft Transmissions 

I— Modelling and Analysis 


Hsiang-Hsi (Edward) Lin 
Memphis State University 
Memphis, Tennessee 

Ronald L. Huston 
University of Cincinnati 
Cincinnati, Ohio 

and 


John J. Coy 
Propulsion Directorate 

U.S. Army Aviation Research and Technology Activity— AVSCOM 
Lewis Research Center 
Cleveland, Ohio 

(NASA-TK-100180) ON DYNAMIC LOADS IN N88-127S7 

FABALLZL SHAF1 IB ANSM ISSICNS . 1: MODELLING 

I AD ANALYSIS (NASA) 16 p CSCL 13l 

unclas 

G3/37 0110985 


December 1987 





US ARMY 
AVIATION 
SYSTEMS COMMAND 
AVIATION R&T ACTIVITY 


E- 3756 


ON DYNAMIC LOADS IN PARALLEL SHAFT TRANSMISSIONS: 


I-MODELING AND ANALYSIS 

Hsiang-Hsi (Edward) Lin* 

Department of Mechanical Engineering 
Memphis State University 
Memphis, Tennessee 38152 

Ronald L. Huston* 

Department of Mechanical and Industrial Engineering 
University of Cincinnati 
Cincinnati, Ohio 45221-0072 

and 

John J. Coy* 

Propulsion Directorate 

U.S. Army Aviation Research and Technology Activity - AVSCOM 

Lewis Research Center 
Cleveland, Ohio 44135-3191 


SUMMARY 

A model of a simple parallel-shaft, spur-gear transmission is presented. 
The model is developed to simulate dynamic loads in power transmissions. Fac- 
tors affecting these loads are identified. Included are shaft stiffness and 
inertia, load and power source inertia, tooth geometry, tooth stiffness, local 
compliance due to contact stress, load sharing, and friction. Governing dif- 
ferential equations are developed and a solution procedure is outlined. A 
parameter study of the solutions is presented in NASA TM-100181 (AVSCOM 
TM-87-C-3) . 


INTRODUCTION 

Recently, there has been increased interest in the dynamic effects in 
gear systems. This interest is stimulated by demands for stronger, higher 
speed, improved-performance, and longer lived systems. This in turn has stimu- 
lated numerous research efforts directed toward understanding gear dynamic phe- 
nomena. However, many aspects of gear dynamics are still not satisfactorily 
understood. 

For example, in industrial settings, a high performance gear system is 
often obtained by overdesigning, and by sacrificing costs, materials, and com- 
pactness. In aerospace and military application where weight is a premium, 
gear systems are often designed under conditions very close to the failure lim- 
its, thereby introducing uncertainties in performance and life prediction. 

They are often prematurely replaced to prevent in-service failure. Moreover, 
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gear systems are often designed by using static analyses. However, when gear 
systems operate at high speed, there are several factors which affect their 
performance. These include shaft torsional stiffness, gear tooth loading and 
deformation, gear tooth spacing and profile errors, rotating speeds, mounting 
alignment, dynamic balance of rotating elements, gear and shaft masses and 
inertia, and the masses and inertias of the driving (power) and driven (load) 
elements . 

There is no agreement among researchers on the best methods for evalu- 
ating dynamic load effects. Hence, gear designers are often confronted with 
conflicting theories. They generally have to rely on past experience, service 
safety factors, and experimental data with a limited range of applicability. 

The objective of this report and NASA TM-1001 81 (AVSCOM TM-87-C-3) is to 
provide more insight into the factors affecting dynamic loads. 

Research efforts on gear system dynamics have been conducted for many 
years. In 1892 Lewis (ref. 1) recognized that the instantaneous tooth load 
was affected by the velocity of the system. In 1925, Earl Buckingham (ref. 2) 
headed an experimental research effort, endorsed by ASME, to measure dynamic 
effects. A report published in 1931 represented the first authoritative docu- 
ment on gear dynamics. It presented a procedure for determining the so-called 
dynamic load increment due to mesh dynamics and gear tooth errors. 

In 1959, Attia (ref. 3) performed an experiment to determine actual 
instantaneous loading. He found that Buckingham's results were conservative. 

In 1958, Niemann and Rettig (ref. 4) found that larger masses caused 
higher dynamic loads, but as the average load became larger the effect of 
larger masses became less important. They also found that very heavily loaded 
gear systems showed no appreciable dynamic load increment, whereas in lightly 
and moderately loaded gear systems there were considerable dynamic load incre- 
ments. In 1958, Harris (ref. 5) suggested that for gear systems isolated from 
external stimuli, there are three internal sources of dynamic loads as follows: 

(1) error in the velocity ratio measured under the working load 

(2) parametric excitation due to stiffness variation of the gear teeth 

(3) nonlinearity of tooth stiffness when contact is lost 

In 1970, Houser and Seireg (ref. 6) developed a generalized dynamic factor for- 
mula for spur and helical gears operating away from system resonances. The 
formula took into consideration the gear geometry and manufacturing parameters 
as well as the dynamic characteristics of the system. 

In 1972, Ichimaru and Hirano (ref. 7) analyzed heavy-loaded spur-gear sys- 
tems with manufacturing errors under different operating conditions. They 
found that the change in tooth profile showed a characteristic trend to 
decrease dynamic load. In 1978, Cornell and klestervelt (ref. 8) presented a 
closed form solution for a dynamic model of spur-gear system and showed that 
tooth profile modification, system inertia and damping, and system critical 
speeds, can have significant effects upon the dynamic loads. In 1981, Kasuba 
and Evans (ref. 9) presented a large scale digitized extended gear modeling 
procedure to analyze spur-gear systems for both static and dynamic conditions. 
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Their results indicated that gear mesh stiffness is probably the key element 
in the analysis of gear train dynamics. They showed that the gears and the 
adjacent drive and load systems can be designed for optimum performance in 
terms of minimum allowable dynamic loads, for a wide range of operating speeds. 

In 1981, Wang and Cheng (ref. 10) developed another dynamic load response 
algorithm. They reported that the dynamic load is highly dependent on the 
operating speed. This model is later modified by Lewicki (ref. 1 1 ) to account 
for the nonlinear Hertzian deformation of meshing gear teeth. The gear dynamic 
load found from the revised model showed little difference from the original 
model since the Hertzian deflection was relatively small in comparison with 
the total gear tooth deflection. Nagaya and Uematsu (ref. 12) stated that 
because the contact point moves along the involute profile, the dynamic 
response should be considered as a function of both the position and speed of 
the moving load. In 1982, Terauchi , et al . , (ref. 13) studied the effect of 
tooth profile modifications on the dynamic load of spur-gear systems. Accord- 
ing to their results, the dynamic load decreased with proper profile 
modi f i cations . 

In this first part of the paper, we present a model of a parallel shaft 
transmission. We consider the effects of shaft stiffness and inertia, load 
and power source inertias, tooth stiffness, local compliance due to contact 
stresses, load sharing, and friction. A parameter study is provided in the 
second part of the paper. 
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NOMENCLATURE 

cross section area of ith element of gear teeth, mm 2 (in. 2) 
damping coefficient, gear tooth mesh, N-sec (lb-sec) 
damping coefficient of shaft, N-m-sec (in. -lb-sec) 
effective modulus of elasticity, N/m 2 (lb/in. 2 ) 
tooth face width, mm (in.) 
shear modulus, N/m 2 (lb/in. 2 ) 

second moment of inertia of ith element of gear teeth, mm 2 (in. 2 ) 

polar moment of inertia of load, m 2 -kg 

polar moment of inertia of motor, m 2 -Kg 

polar moment of inertia of gear 1, m 2 -Kg 

polar moment of inertia of gear 2, m 2 -Kg 

stiffness of gear tooth, N-m/rad 

stiffness of shaft, N-m/rad 
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distance between elements i and j , mm 

gear tooth deformation due to beam deflection, mm 

gear tooth deformation due to contact deformation, mm 

gear tooth deformation due to foundation flexibility, mm 

total gear tooth deformation, mm (in.) 

base radii of gears, mm (in.) 

pitch radi i of gears , mm 

thickness of element i, mm 

torque on load, N-m 

torque on motor, N-m 

torque on gear 1 , N-m 

torque on gear 2, N-m 

sliding velocity during tooth mesh, mm/sec 

rolling velocity, mm/sec 

applied load, N/m (lb/in.) 

x-coordinate of element i, mm (in.) 

y-coordinate of element i, mm (in.) 

load angle, rad 

backlash, mm (in.) 

angular displacement, rad 

angular velocity, rad/sec 

angular acceleration, rad/sec 2 

lubricant viscosity, N-sec/m 2 (lb-sec/in. 2 ) 

Poisson's ratio 
damping ratio 


MODELING 

Figure 1 depicts a model of the transmission. It consists of a motor or 
power source connected by a flexible shaft to the gear system. The gear sys- 
tem consists of a pair of involute spur gears. They are connected to the load 
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by a second flexible shaft as shown. Symbolically, the model may be repre- 
sented by a collection of masses, springs, and dampers as in figure 2. 

Let ©m, ©i, © 2 , and ©|_ represent the rotations of the motor, the 
gears, and the load. Then by using standard procedures of analysis, the gov- 


erning differential equations for the rotations may be written as 

Jm§m + C s i (©m - ©l ) + K s i (©m - ©i > = Tm (1) 

J 1 ©1 + C s i<©i - ©m> + i <©i - ©m> + Cg(t)[Rb]©i - Rb2©2^ 

+ Kg(t)[Rbl (Rbl®l - Rb2©2^ = Tf i C t ) (2) 

J2©2 + C S 2(©2 - ©1 > + Ks 2^®2 - ©l ) + Cg(t)[Rb2©2 ~ Rbl©ll 

+ Kg(t)CRb2(Rb2©2 - Rbl©l>3 - T f2 <t> (3) 

J|_©L + Cs2^©L ~ ©2^ + ^s2^©L ~ ©2^ = -T|_ (4) 


where Ji , J2» ancl J L represent the mass moments of inertia of the 
motor, the gears, and the load; C S ] , C S 2> and Cg(t) are damping coefficients 
of the shafts and the gears; K s -j , K S 2 , and K g (t) are stiffnesses of the 
shafts and the gears; T^, Tl, Tfi(t), and Tf 2 (t) are motor and load torques 
and frictional torques on the gears; Rbi and Rb2 are base circle radii of 
the gears; t is time; and the dots over © indicate time differentiation. 

In developing equations (1) to (4) several simplifying assumptions are 
employed: 

(1) The dynamic process is studied in the rotating plane of the gears. 
Out-of-plane twisting and misalignment are neglected. 

(2) Damping due to lubrication of the gears and shafts is expressed in 
terms of constant damping factors. 

(3) The differential equations of motion are developed by using the theo- 
retical line of action. 

(4) Low contact ratio gears are used in the analysis. Specifically, the 
contact ratio is taken between 1 and 2. 


ANALYSIS 

A major task in the analysis is to determine the values of the stiffness, 
damping, and friction coefficients appearing in equations (1) to (4). Another 
task is to determine the ratio of load sharing between the teeth during a mesh 
cycle. These factors depend on the roll angle of the gears. Thus, equations 
(1) to (4) are made nonlinear by these terms. 
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Stiffness 


Gear stiffness . - Consider first the stiffness coefficients kg, K S ] , 
and k S £. Let the tooth surface have the form of an involute curve. Let Wj 
be the transmitted load at a typical point j of the tooth profile. Let qj 
be the deformation of the tooth at point j in the direction of W j . Then 
the gear stiffness kg j for the gear teeth in contact at j is 




(5) 


In general qj will depend on the following: (1) the bending of the 

tooth, the shear deformation of the tooth, and the axial compression of the 
tooth; (2) the deflection due to the flexibility of the tooth foundation and 
fillet; and (3) the local compliance due to the contact stresses. 

To determine qj let the tooth be divided into elements as shown in 
figure 3. Let i be a typical element with thickness T ^ , cross section 
area , and second moment of inertia Ii- Let L^j be the distance between 
element i and point j along the x-axis. Let |3j be the angle between 
Wj and the y-axis. (See fig. 3.) 

Consider the tooth to be a nonuniform cantilever beam. Let qpj be the 
contribution to qj by the bending, shear, and axial deformation of the 
tooth. Then q^j may be represented as the sum of the deformation in the ele- 
ments i beneath point j. That is. 


n 



i-1 


( 6 ) 


where n is the number of elements beneath j and qbi j is the deformation 
of element i due to the load Wj. 


With standard analysis qbn is found to have the value (refs. 8, 14, 
and 15): 


"mi * ( W 


2 

cos |3 . 
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- cos |3. sin |3. 
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j ♦ 2 Vj L u)/ 2I i] 

cos 2 |3j[l2(l + v)T./5A.] + sin 2 |3j(T. /A. ) J 


( 7 ) 


where is the half-tooth thickness at element i (see fig. 3), v is 

Poisson's ratio, and E e is the "effective elastic modulus" depending upon 
whether the tooth is wide (plane strain) or narrow (plane stress). Specifi- 
cally, for a "wide" tooth, where the ratio of the width to thickness at the 
pitch point exceeds 5 (ref. 14), E e is 
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E 


( 8 ) 



where E is Young's modulus of elasticity. For a "narrow" tooth (width-to- 
thickness ratio less than 5), E e is 


E e = E (9) 

Expressions similar to equations (7) hold for qf j the contribution to 
qj for the deformation due to the flexibility of the tooth fillet and founda- 
tion (ref. 15). 


Let q C j be the contribution to qj from the local compliance due to 
contact stresses. With the procedures of Lundberg and Palmgren (ref. 16) q c j 
may be expressed as 



1 .275 

E 0 . 9 F 0 . 8 W 0.1 


(10) 


where F is the width of the tooth. 

Hence, by superposition, the deformation at j in the direction of Wj 


is 


Qj - q bj ♦ Q fj ♦ q cj 

The above expressions were used to calculate the deformations for two dif- 
ferent gear pairs. The results are shown graphically in figures 4(a) and (b). 


Shaft stiffness . - The shaft stiffness K s is given by the standard 
expression 


K 



( 12 ) 


where G is the shear modulus, i is the shaft length, and 3 is the polar 
moment of area given by 


J 



(13) 


where D i s the shaft diameter. 


Damping 

Shaft damping . - Next consider the damping coefficients C s i, C S 2 , and 
Cg- Damping in the shafts is due to the shaft material. In equations (1) to 
(4) the coefficients C s i and C S 2 are taken to have the form 
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/ 


( 14 ) 


(15) 


where 5 S represents the damping ratio. Experiments have shown that £ s has 
values between 0.005 and 0.075 (ref. 17). 


Mesh damping . - Similarly, the effect of damping of the gear mesh is 
taken as 


1/2 


C„ = 2S 


R bl J l * R b2 J 2 


(16) 


where, as before, % is the damping ratio. Measurements have shown £ to 
have values between 0.03 and 0.17 (refs. 9 and 10). 

Friction . - Equations (1) to (4) contain terms Tfi and Tf 2 which rep- 
resent the frictional moments of the driving and driven gears. These moments 
occur because of the relative sliding of the gear teeth. Buckingham (ref. 18) 
has recorded a semiempirical formula for the friction coefficient f of 
boundary lubrication as 


f = 0.05e‘°- 125Vsa + 0.002 (17) 

where V S £ is the sliding speed measured in in. /sec. An analogous expression 
for elastohydrodynamic lubrication has been developed by Benedict and Kelly 
(ref. 19) and by Anderson and Loewenthal (ref. 20) as follows: 

f = 0.0127 log f45.94 W/F^V^vO (18) 


where 

W the applied load, N/m (Ib/in.) 

F face width, mm (in. ) 

Vr rolling velocity, mm/sec (in. /sec) 
y 0 lubricant viscosity, N-sec/m 2 (lb-sec/in. 2) 
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Figures 5(a) and (b) show graphs of the friction coefficient as given by 
equations (17) and (18) as a function of the roll angle. Figures 6(a) and (b) 
show the resulting effect upon the friction torque. 


Mesh Analysis 

Figure 7 illustrates the motion of a pair of meshing teeth. The initial 
contact occurs at A, where the addendum circle of the driven gear intersects 
the line of action. As the gears rotate the point of contact will move along 
the line of action APD. When the tooth pair reaches B, the recessing tooth 
pair disengages at 0 leaving only one zone. When the tooth pair reaches 
point C, the next tooth pair begins engagement at A and starts another 
cycle. 

In the analysis, the position of the contact point of the gear teeth 
along the line of action is expressed in terms of roll angles of the driving 
gear tooth. 

Figure 8 shows typical stiffness and load sharing characteristics through 
a mesh cycle. Let a series of mating tooth pairs be denoted as a, b, c, d 
and let points A, B, P, C, D be the same as those in figure 7. Then, AB 
and CD represent the double contact regions, BC represents the single con- 
tact region, and as before P is the pitch point. 

The stiffness values at double contact regions are clearly much higher 
than those at single contact regions. When gears rotate at appreciable speed, 
this time-varying stiffness as shown in figure 8 is the major excitation 
source for the dynamic response of the system. 


DISCUSSION 

The objective of this analysis is to establish the governing differential 
equations and to present a procedure for solution. As noted, the equations 
themselves are nonlinear. However, they may be efficiently solved by using a 
1 inearized-iterati ve procedure as follows: 

The linearized equations may be obtained by dividing the mesh cycle into 
n equal intervals. Let a constant input torque T^ be assumed. Let the out- 
put torque T[_ be fluctuating because of damping in the gear mesh, because of 
friction, and because of time-varying mesh stiffness. 

Let initial values of the angular displacements be obtained by preloading 
the input shaft with the nominal torque carried by the system. Initial values 
of the angular speeds may be taken from the nominal operating speed of the 
system. 

The iterative process is then as follows: the calculated values of the 

angular displacements and angular speeds after one period are compared with 
the assumed initial values. Unless the differences between them are suffi- 
ciently small, the procedure is repeated by using the average of the initial 
and calculated values as new initial values. 
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Finally, observe that the term < Rb 1 - R b2 e 2^ i n the equation of motion 
represents the relative dynamic displacement of the gears. Let S represent 
the backlash. Let gear 1 be the driving gear. The following conditions can 
occur 


Case 1 


The normal operating case is 

R bl e l - R b2 e 2 > 0 

The dynamic mesh force F is then 

F = Kg(t)(Rbi©-| - Rb2©2^ + Cg < t ) < R^ ] ©i - Rb2®2^ 


(19) 

( 20 ) 


Case 2 

R bl e l " R b2 e 2 = 0 and | Rbl e l ~ R b2 e 2 I = & (21) 

In this case, the gears will separate and the contact between the gears will 
be lost. Hence, 

F = 0 (22) 


Case 3 

R bl 9 l - R b2 e 2 < 0 and I R bl ®1 - R b2 e 2l > $ (23) 

In this case, gear 2 will collide with gear 1 on the backside. Then, 

F = Kg(t)[(Rb2©2 ~ R bl®P - S] + Cg(t)(Rbi©i - Rb2®2^ (24) 


CONCLUSION 

A low contact ratio spur-gear transmission model is developed. The model 
includes inertias of load and power source, stiffness of shaft, time-varying 
mesh stiffness, and damping and friction inside gear transmissions. 

Governing equations of the model are derived and a linearized iterative 
procedure for solution is presented. Parameter study including rotating 
speed, diametral pitch, applied load, damping, stiffness, and inertia will be 
presented in NASA TM-1001 81 (AVSCOM TM-87-C-3). 
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FIGURE 1. - A SIMPLE SPUR GEAR SYSTEM. 


FIGURE 2. - A MATHEMATICAL MODEL OF THE TRANSMISSION SYSTEM. 



FIGURE 3. - ELEMENT MODELLING OF A GEAR TOOTH. 
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NORMALIZED DEFORMATION, qEF/W 




ROLL ANGLE, DEG 

(B) PITCH RADII, R p1 = 57.1 MM, R p2 - 171.3 MM. 

FIGURE 4. - NORMALIZED DEFORMATION OF A PAIR OF 
TEETH. MODULE, 3.18 MM; PRESSURE ANGLE, 20°; 
FACE WIDTH, F = 25. *1 MM; MODULUS OF ELASTICITY, 
E = 207 GPa; APPLIED LOAD, W = 105 kN/m. 



(B) BENEDICT AND KELLEY'S MODEL. 

FIGURE 5. - FRICTION COEFFICIENT OF GEARS IN FIG- 
URE 4(A) AT 1500 RPM. 
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TORQUE, N-M 




(B) BENEDICT AND KELLEY'S MODEL. 

FIGURE 6. - FRICTION TORQUE VARIATION ALONG THE CON- 
TACT PATH FOR GEARS IN FIGURE 4(A) AT 1500 RPM. 



FIGURE 7. - ILLUSTRATION OF GEAR MESHING ACTION. 




0 1 2 3 4 5 

NORMALIZED CONTACT POSITION (BASE PITCH) 


FIGURE 8. - TYPICAL STIFFNESS AND LOAD SHARING OF 
LOW CONTACT RATIO GEARS OF FIGURE 4(A). 
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